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m 要 : 土壤 盐 结 皮 对 干旱 区 土壤 演化 和 生态 水 文 过程 具 有 重要 影响 ,目前 的 研究 对 于 盐 结 皮 土 壤 
的 水 盐分 布 特征 涉及 较 少 ,并 且 未 考虑 盐 结 皮 厚 度 的 影响 ,导致 研究 结果 存在 较 大 差异 。 因 此 , 通 
过 室内 模拟 试验 ,设置 4 种 初始 盐分 浓度 (0 g L7 ,10 g* L',150 g- L $28 250 g- L) DLZR SUR [8] 35 Zi 
厚度 (4.5 mm、6.6 mm ff 7.3 mm), RA 38 22 38 4 COP 3. EN 3E LALA IC B 20 CE 
tho BRERA) 对 比 无 盐 处 理 , 含 盐 处 理 盐 结 皮 越 厚 , 土 壤 剖 面 含水 量 越 大 , 伪 盐 量变 化 幅度 越 


小 。(2) 试验 结束 时 ,4.5 mm 厚度 盐 结 皮 土 壤 含 水 量 


度 盐 结 皮 土壤 剖面 含水 量 显 著 大 于 无 盐 处 理 (P<0. 


分 布 特征 和 无 盐 处 理 相似 ,6.6mm 和 7.3 mm 厚 
05)。(3) 试验 结束 时 4.5 mm 6.6 mm 和 7.3 mm 


厚度 盐 结 皮 土 壤 训 面 最 小 含 盐 量 相 比 初 始 含 盐 量 分 别 减 小 了 90.5% 46.3% 和 32.1%。 研 究 结 果 验 


~ 


相关 研究 需 综 合 考虑 盐 结 皮 厚度 的 影响 。 
关键 词 : 盐 结 皮 厚 度 ; 水 盐分 布 特征 ; RAE; HYDRU 
文章 编号 : 1000 - 6060(2023)08 - 1303 - 11(1303 ~ 1313) 


盐 结 皮 是 土壤 盐 省 化 的 极端 表现 ,对 土壤 性 质 
以 及 地 表 水文 过 程 等 具有 重大 影响 FAKE 
壤 水 盐 运 动 是 一 个 复杂 的 动态 过 程 ”, 并 且 土 壤 水 
盐分 布 特征 会 影响 土壤 物理 化 学 及 生物 过 程 “ , 进 
影响 土壤 重金 属 再 分 布 “\ 碳 排放 "等 ,而 在 干旱 
区 盐 结 皮 普 遍 存 在 ,因此 人 研究 盐 结 皮 土 壤 水 盐分 布 
特征 具有 重要 的 理论 意义 。 
目前 国内 外 学 者 对 于 土壤 水 盐分 布 的 研究 主 
要 集中 在 非 盐 结 皮 土壤 。Zhang 等 5 对 蒸发 过 程 中 
浅 层 裸 土 水 盐分 布 进 行 室内 试验 ,发 现 由 于 水 盐 迁 
移 不 同步 ,盐分 剖面 呈 波 浪 形 分 布 。 当 土壤 表面 被 
不 同 材料 履 盖 时 会 显著 降低 土壤 蒸发 ,减轻 土壤 盐 
分 的 表 聚 ”%"。 盐 结 皮 对 土壤 的 影响 类 似 于 秸秆 覆 
盖 ” ,因此 盐 结 皮 必然 会 影响 土壤 水 盐分 布 特征 。 
然而 ,由 于 目前 的 技术 手段 对 高 盐 区 域 土壤 含 盐 量 
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证 了 擂 结 皮 厚 度 会 对 土壤 水 盐分 布 产生 较 大 影响 ,因此 建议 未 来 对 于 盐 结 皮 土 壤 水 盐分 布 特征 的 


S-1D 模型 


测量 的 局 限 性 以 及 盐 结 皮影 响 下 土壤 水 盐 热 耦合 
作用 的 复杂 性 ,导致 目前 仅 有 少数 学 者 对 盐 绪 皮 土 
壤 水 盐分 布 特征 进行 研究 ,如 Gran 等 ”通过 室内 试 
验 ,发 现 盐 结 皮影 响 下 不 同 质地 土壤 含水 量变 化 曲 
线 有 很 大 差异 ,盐分 在 0~4 cm 土 层 累积 ,并 且 越 接 
近 土 壤 表 面 盐 分 浓度 越 高 。 英 治 新 等 “研究 表明 
盐 结 皮 土 壤 含 水 量 均 大 于 无 盐 结 皮 土 壤 , 不 同 盐 结 
皮 厚 度 会 影响 水 分 的 迁移 导致 不 同 土 层 含水 量变 
化 规律 不 同 ,并 且 盐 结 皮 越 厚 的 土壤 区 域 其 各 层 土 
寺中 的 盐分 含量 越 高 。 王 晓 静 等 "对 塔里木 沙漠 
公路 防护 林 市 土壤 进行 分 析 , 发 现 盐 绪 皮 厚 度 会 逐 
渐 增 长 变 厚 并 趋 于 稳定 ,土壤 电导 率 呈 现 先 增加 后 
逐渐 降低 的 趋势 。 Guglielmini 等 56 Dashtian $^" ah 
过 建立 模型 预测 盐分 在 土壤 中 的 迁移 过 程 ,结果 表 
明 ,盐分 浓度 最 大 值 出 现在 液 - 汽 界面 附近 。 这 一 
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液 - 汽 转化 面 被 称 为 蒸发 锋 “” ,能 直接 反映 土壤 
水 盐分 布 状 况 ,其 变化 会 影响 盐分 累积 的 深度 和 水 
分 运动 状态 ,是 描述 蒸发 过 程 中 水 盐分 布 的 重要 指 
标 ”。 然 而 ,当前 研究 针对 盐 结 皮 土 壤 蒸 发 锋 提出 
的 不 同 观点 ,使 盐 结 皮 土 壤 水 盐分 布 特征 的 研究 仍 
然 存在 争议 。 一 些 学 者 ” 认为 ,蒸发 过 程 中 盐 结 
土壤 蒸发 锋 区 域 较 窗 (2~4 cm) ,而 男 一 些 学 者 “™ 
在 研究 中 发 现 有 较 宽 的 蒸发 锋 区 域 (8~12 cm)。 上 
述 研 究 结 果 的 差异 ,可 能 原因 是 不 同 的 试验 条 件 ， 
如 :盐分 浓度 ””、 盐 分 类 型 ”不 同 地 下 水 深度 
等 ,导致 形成 不 同 厚度 的 盐 结 皮 , 从 而 对 土壤 水 盐 
分 布 特征 产生 不 同 的 影响 ,但 上 述 研 究 均 未 考虑 盐 
结 皮 厚度 这 一 重要 的 物理 特征 。 

在 干旱 半 干 旱地 区 土壤 盐 溃 化 程度 不 同 , 其 土 
壤 含 盐 量 有 很 大 的 空间 变异 性 ,各 种 厚度 盐 结 皮 广 
沁 分 布 ,而 盐分 浓度 是 影响 盐 结 皮 形 成 发 育 的 重要 
因素 ””。 人 研究 表明 ,蒸发 是 土壤 水 盐 运 动 的 主要 
驱动 力 ” ,而 不 同 盐 结 皮 厚 度 会 使 土壤 蒸发 产生 较 
大 差异 *” ,因此 盐 结 皮 的 厚度 必然 会 影响 土壤 水 
盐分 布 特征 ,但 目前 盐 结 皮 厚 度 对 土壤 水 盐分 布 特 
征 的 作用 机 理 还 不 明确 。 因 此 ,本 人 研究 设置 了 4 种 
不 同 初始 浓度 以 得 到 不 同 厚度 的 盐 结 皮 ,分析 其 土 
壤 水 盐 的 动态 变化 ,明确 不 同 厚 度 盐 结 皮 土 壤 水 盐 
分 布 特征 的 差异 ,阐明 盐 结 皮 厚 度 对 土壤 水 盐分 布 
特征 的 影响 机 理 , 为 定量 描述 干旱 区 盐 涡 土壤 的 地 
表 水 文 过 程 提供 理论 基础 。 


1 材料 与 方法 


1.1 试验 材料 

供 试 土壤 质地 为 砂 士 ,土壤 容重 为 1.65 e-em, 
饱和 含水 率 为 39.67%。 采 用 直径 10 em 33 em 的 
HR ls CPVCÓE EIN Gus ,并 包 于 隔 热 材料 。 将 
试验 采用 的 砂 土 洗 净 ,去 除 杂 质 并 将 电导 率 控制 在 
100 uS- cm- 以 下 , 烘 干 后 过 0.5 mm + Be 85535] 5] Re 
入 PVC 管 ,在 PVC 管 底部 填充 3 em 高 的 滤 石 (直径 
5 mm)。 试 验 土 柱 自 下 而 上 充分 饱和 24 h。 在 土 柱 
正 上 方 38 em IIb JJCEL— se d RAT (50 W) H PIKSI 
发 ,整个 试验 过 程 平均 环境 温度 为 25+3.5 ,湿度 
为 (47.2+13.9)%。 
1.2 试验 设计 与 方法 

试验 设置 了 4 个 人 处理, 分 别 为 0 g:L'(CK)、10 


g-L'(LC),150 g-L'(MC)f250 g: L' (HC), EP 
CK 为 无 盐 处 理 , 采 用 蒸馏 水 饱和 ,LC、MC 和 HC 采 
用 不 同 浓度 NaCl 溶 液 进行 初始 饱和 ,以 获得 不 同 厚 
度 的 盐 结 皮 。 试 验 蒸 发 采用 自动 称 重 平 台 ( 型 号 
WP20 ,北京 时 域 通 科技 有 限 公 司 ) 测 定 ,并 连接 数据 
采集 髓 (型 号 i-logger, 北 京 时 域 通 科 技 有 限 公 司 ) 记 
录 数 据 。 采 取 部 分 重复 逐步 退出 的 方法 3 , 即 在 选 
定 的 时 间 点 将 部 分 土 柱 进行 切割 ,将 盐 结 皮 从 土壤 
表面 剥离 ,10 em 深度 以 上 每 隔 2 cem 取 土 ,10~30 cm 
每 隔 $ cm 取 土 ,对 其 含水 量 及 含 盐 量 进行 测定 。 根 
据 之 前 研究 关 所 得 不 同 初始 浓度 各 蒸发 阶段 的 时 
间 ,本 文选 取 5 个 关键 节点 作为 退出 时 间 ,其 中 第 一 
阶段 选取 特征 时 间 a, 第 二 阶段 选取 特征 时 间 b 和 c， 
第 三 阶段 选取 特征 时 间 d 和 e, 不 同 处 理 各 特征 时 间 
如 图 1 所 示 。 采 用 烘 干 法 测定 土壤 含水 量 ,将 烘 干 
土 样 赋 及 .过 1 mm 土壤 得 , 按 填 水 比 为 1 配置 成 深 
液 , 用 玻璃 棒 搅 拌 . 静 置 ,之 后 使 用 电导 率 仪 测定 其 
电导 率 , 计 算 土壤 含 盐 量 与 电导 率 之 间 的 关系 式 如 
公式 (1) 所 示 。 采 用 游标 卡尺 测量 盐 结 皮 厚度 ,之 
后 将 盐 结 皮 加 蒸馏 水 溶解 .过滤 、 烘 干 后 测量 盐 结 
皮质 量 。 


S, =4.533 x EC - 0.723 (1) 
AP: S HERSE (eke!) SEC A IK EUG 1:5 
提 液 电导 率 (mS.cm-)。 
土壤 进 气 值 指 饱 和 土壤 脱水 过 程 中 开始 进入 
空气 时 的 吸力 值 ,与 Van Genuchten 模 型 中 进 气 吸力 
相关 参数 和 形状 系数 有 关 “ 。 本 研究 通过 引用 土 
壤 进 气 值 来 分 析 不 同 处 理 的 水 力学 参数 。 土 壤 进 
气 值 的 计算 公式 如 下 : 


wd) Me Cu (2) 


式 中 :万 为 土壤 进 气 值 ; a n 分 别 为 VG 模型 中 与 进 
气 值 有 关 的 经 验 参数 。 

受 试验 条 件 限制 ,所 有 人 处理 蒸发 锋 的 位 置 变 化 
及 下 移 时 间 并 不 能 被 完全 观测 ,因此 选择 HYDRUS- 
1D 模 型 延长 蒸发 时 间 ,继续 观测 土壤 水 分 动态 过 
程 ,以 明确 各 处 理 蒸发 锋 的 变化 。 
13 土壤 水 盐 动 态 模型 

HYDRUS-1D 模 型 中 土壤 水 分 运动 方程 : 


36(h) a(h) 
3-7 exe (n + J (3) 
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式 中 :0(h) 为 土壤 体积 含水 率 (cm.cm”);h 为 压力 水 
$ (em) ;为 水 分 运 移 时 间 (d) ;z 为 垂直 方 问 坐标 变 
tg Cem) ; K(h) NARA ZK A (em - d! ) 


盐分 运 移 方程 : 
AC) _ af anac |. Age) 
at = iic dz s 


式 中 :C 为 浴 质 浓度 (gecm );D 为 弥散 系数 (cm*…d'); 
z; HE ETT [Hy At (em) 5 q 297 GE E (em hi!) 
模型 高 度 为 30 em, 412 31 P 75 8, 73 5 JR] ER 
1 cm ,分 别 在 土 层 2 cm .4 em,6 cm 、 8 cm 10 em, 15 
cm、20 cem、25 cm 和 30 cm 处 设 观测 点 。 模 型 模拟 时 
间 同 试验 时 间 相 同 ,之 后 延长 蒸发 时 间 到 120 d , 49] 


^ 特征 时 间 


ZE RE mm:d' 


时 间 /d 


RRS md 


时 间 /d 


始 时 间 步 长 为 0.01 d, 最 小 时 间 步 长 为 0.001 d, 最 大 
时 间 步 长 为 1 ds 

通过 HYPROP 软件 实测 土壤 水 分 特征 曲线 得 
到 各 处 理 土壤 水 力 参 数 ,溶质 运 移 参数 基于 文献 确 
定 2 2 ,之 后 通过 HYDRUS-1D 模型 采用 实测 土壤 
含水 量 与 含 盐 量 对 参数 进行 修正 ,最 终 不 同 处 理 土 
壤 水 力 参 数 与 溶质 运 移 参数 如 表 1 所 示 。 本 试验 中 
采用 同 种 土壤 质地 ,不 同 初始 盐分 浓度 饱和 土壤 ， 
导致 各 处 理 土壤 含 盐 量 有 很 大 不 同 ,而 土壤 盐分 含 
量 会 影响 土壤 持 水 性 和 土壤 孔隙 ,从 而 使 各 处 理 
土壤 水 力 参 数 产 生 较 大 差异 , 随 着 初始 浓度 的 增 
加 ,a 逐渐 减 小 ,n 逐 渐 增 大 1。 


蒸发 曲线 


RATA md | 


时 间 /d 


蒸发 速率 /mm"d7 


时 间 /d 


注 :CK 为 0g:L';LC 为 10 g-L'; MC H 150 gL’; HC 29250 gL';a~e 为 试验 特征 时 间 。 下 同 。 
图 1 各 处 理 逐 步 退出 时 间 节 点 的 蒸发 速率 


Fig. 1 Evaporation rate of each treatment at gradually exit time nodes 


表 1 不 同 处 理 土壤 水 力 人 参数 


Tab.1 Soil hydraulic parameters under different treatments 


不 同 浓度 人 处理/g*1L" 0, /cm cm 0. /cm* cm? 
0 0.0029 0.420 
10 0.0029 0.396 
150 0.0030 0.396 
250 0.0030 0.396 


a n k/cm*d" l 
0.0320 3.80 440.67 0.5 
0.0178 6.60 440.67 0.5 
0.0057 19.00 440.67 0.5 
0.0053 19.45 440.67 0.5 
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1.4 模型 初始 条 件 和 边界 条 件 

对 于 土壤 水 分 运动 ,试验 开始 时 土壤 初始 饱 
和 ,初始 条 件 为 饱和 含水 量 , 上 边界 设置 为 变通 量 
边界 ,将 试验 实测 的 蒸发 通 量 作 为 上 边界 ,下 边界 
为 零 通 量 边界 。 对 于 溶质 运 移 ,不 同 处 理 初始 条 件 
为 初始 饱和 浓度 ,上 下 边界 条 件 均 设 置 为 浓度 通 量 
边界 。 
1.5 数据 处 理 

采用 Excel 2019 Fl SPSS 25 软件 进行 统计 分 析 ， 
使 用 Origin 2018 软件 绘制 图 形 。 模 型 验证 使 用 决 
ERAUR MEINIR (RMSE) , RMSE 越 接近 于 
0, RRR AIER ER fey, RC pem 1, WERA 
以 较 好 地 捕捉 到 实测 值 的 变动 趋势 。 


(5) 
式 中 :NN 为 样本 容量 ;5; 为 模拟 值 ;M; 为 实测 值 。 


2 结果 与 分 析 


2.1 不 同 厚 度 盐 结 皮影 响 下 土壤 水 分 动态 变 

图 2 所 示 为 试验 结束 时 含 盐 处 理 土壤 盐 结 
皮 。 含 盐 处 理 在 不 同 特征 时 间 盐 结 皮 的 质量 与 厚 
度 存 在 差异 (图 3) , 含 盐 处 理 盐 结 皮 随 着 时 间 的 
推移 质量 逐渐 增加 ,厚度 逐渐 变 厚 ,土壤 含 盐 量 对 
盐 结 皮 厚 度 的 影响 较 大 ” 。 李 胜 辉 等 ”研究 发 
现 盐 浓度 从 0.5 mol + L 34705] 2.5 mol + LAY , 盐 结 
皮 厚 度 逐 渐 增 加 ,本 试验 研究 结果 与 其 相似 ,在 试 
验 结 束 时 LC、MC 和 HC 盐 结 皮 厚 度 分 别 为 4.5 mm, 
6.6 mm 和 7.3 mm。 通 过 分 析 不 同 处 理 土壤 含水 量 
随时 间 的 动态 变化 (图 4) , MCA HC 土壤 含水 量 相 


(a) LC 


N 


比 CK 有 显著 差异 (P<0.05)。 对 于 本 研究 定义 的 表 
层 土壤 (0~2 em) , 盐 结 皮 越 厚 土 壤 含 水 量 越 大 。 在 
特征 时 间 a,CK、LC、MC A HC 土壤 含水 量 分 别 为 
0.123 cm? .cm 0.202 cm? .cm ?,0.329 cm .cm 和 
0.293 cm em ,这 是 因为 土壤 中 盐 溶 液 使 士 壤 渗透 
势 发 生变 化 ,饱和 水 汽 压 减 小 ” ,抑制 水 分 的 向 上 传 
输 ,导致 仿 盐 处 理 含水 量 大 于 CK。 在 特征 时 间 b， 
HC 土壤 含水 量 显 闭 性 高 于 CKCP<0.05) , 相 比 特征 
时 间 a 仅 减少 了 0.019 cm em? ,可 能 原因 是 此 时 HC 
土壤 表面 形成 一 层 较 薄 的 盐 结 皮 , 抑 制 了 水 分 同上 
传输 ,导致 土壤 含水 量变 化 较 小 。 在 特征 时 间 e,CK 
表层 土壤 含水 量 降 低 到 0.009 em? - em? ,表层 土壤 基 
本 干燥 (土壤 合 水 量 小 于 风干 土 含水 量 0.050 cm- 
cm 定义 为 干燥 土 ”) ,而 含 盐 处 理由 于 盐 结 皮 逐 渐 
变 厚 导致 水 分 传输 阻力 的 增加 ,此 时 LC、MC 和 HC 
土壤 含水 量 分 别 为 0.096 cm?+cm™ ,0.194 en? cem? fll 
0.261 cm cm“”。 在 特征 时 间 d~e 阶段 ,各 人 处理 表层 
土壤 含水 量 仅 有 小 幅 变 化 ,试验 结束 时 LC、MC 和 
HC 土壤 舍 水 量 均 大 于 CK(0.005 em? em?) ,分 别 为 
0.008 cm^* cm ^,0.131 cm cm 和 0.176 cm? em ^, ER 
结 皮 越 厚 ,表层 十 壤 含 水 量 越 大 。 

对 于 本 研究 定义 的 上 层 土壤 (2~10 em) , 含 盐 
处 理 土 壤 含水 量 在 该 层 变 化 幅度 较 大 。 在 特征 时 
HJ a~c 阶段 ,由 于 表层 水 分 的 减少 ,上 层 土壤 水 分 开 
始 向 表层 迁移 ,HC 由 于 较 低 的 蒸发 驱动 力 ”“, 叶 致 
其 上 层 土壤 含水 量 相 比 其 他 处 理 变化 较 小 。 在 特 
征 时 间 d~e 阶 段 ,LC 其 2~4 em 土 层 含水 量 相 比 0~ 
2 cm 突然 增加 ,土壤 含水 量 的 突变 表明 蒸发 锋 移 动 
到 了 2~4 cm 深度 , 当 蒸 发 锋 发 生 明 显 下 移 , 其 上 方 
土壤 水 力 联 系 被 切断 ,此 时 LC 土壤 表面 盐 结 皮 厚 度 


(b) MC (c) HC 


图 2 试验 结束 时 含 盐 处 理 土壤 盐 结 


Fig.2 Soil salt crust on saline treated at the end of the test 
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图 3 不 同 初 始 浓度 盐 结 皮 质量 与 厚度 
Fig.3 Mass and thickness of salt crust with different initial 


concentrations 
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m 19 INC [—wcL HC x 与 质量 趋 于 稳定 ,这 与 Nachshon 45 ^ Br ZH RAB 
aE FF H SW. MMC 和 HC 上 层 土壤 含水 量 均 大 于 0.14 
30 ae 1 = em' cm , 仍 与 表层 土壤 保持 水 力 联系 ,未 明确 观测 
"E B i = 到 蒸发 锋 的 下 移 , 主 要 原因 是 试验 前 期 土壤 洲 液 中 
E 20 Lc xinh HC BES AER SF xU EK SP IS RAS IF 
i - = 盐 结 皮 的 快速 形成 阻碍 了 水 分 的 传输 ,导致 土壤 含 
各 m 水 量 相 比 CK 变 化 较 小 。 试 验 结 束 时 ,上 层 土壤 含 


水 量 表现 为 HC>MC>LC>CK, 盐 结 皮 越 厚 , 上 层 土壤 
含水 量 越 大 。 

对 于 下 层 土壤 (10~30 em) ,不 同 处 理 均 表现 出 
随 着 深度 的 增加 土壤 含水 量 逐 渐 增 加 。 在 试验 周 
期 内 , 含 盐 处 理 下 层 土 壤 与 CK 相 比 变化 较为 平缓 ， 
盐 结 皮 越 厚 , 下 层 土壤 含水 量 越 大 ,主要 原因 是 随 
着 时 间 的 推移 , 盐 结 皮 逐 渐 生 长 变 厚 ,导致 蒸发 驱 


土壤 含水 量 /cms'cm- 
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图 4 不 同 处 理 各 深度 土壤 含水 量 随时 间 的 变化 


Fig. 4 Variation of soil water content at different depths with time in different treatments 
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图 5 不 同 处 理 各 深度 土壤 含 盐 量 随 时 间 的 变化 


Fig. 5 Variation of soil salt content at different depths with time in different treatments 


动力 降低 ,下 层 土壤 水 分 向 上 迁移 速率 减缓 。 在 试 
验 结束 时 ,LC 下 层 土壤 含水 量 相 比 CK 较 小 ,而 MC 和 
HC 显著 大 于 CK(P<0.05), 可 能 原因 是 在 试验 后 期 
LC 盐 结 皮 较 薄 , 对 土壤 水 分 的 抑制 作用 远 小 于 MC 和 
HC, Jf ERA BOE jS TAR I TRU ,在 一 
定 程度 上 增加 了 蒸发 速率 ” ,这 一 特殊 效应 导致 LC 
下 层 土壤 含水 量 小 于 CK ,而 对 于 HC 和 MC ,这 一 特 
殊 效 应 未 对 盐 结 皮 的 抑制 作用 产生 明显 影响 ,从 而 
导致 HC 和 MC 下 层 土壤 含水 量 远 大 于 CK。 
22 不 同 厚度 盐 结 皮影 响 下 土壤 盐分 动态 变化 

图 $ 所 示 为 含 盐 处 理 各 深度 土壤 含 盐 量 随时 间 
的 变化 。 对 于 表层 土壤 ,不 同 处 理 盐分 均 在 土壤 表 
面 累积 ,这 是 因为 土壤 中 盐 涂 液 在 蒸发 驱动 下 逐渐 
向 土壤 表面 迁移 , 随 着 溶液 中 水 分 的 蒸发 ERO 
留 在 土壤 表面 , 当 盐 分 浓度 达到 饱和 状态 时 ,盐分 开 
始 在 土壤 表面 结 品 析出 。 在 特征 时 间 a,LC 表 层 土壤 
含 盐 量 大 于 初始 值 (3.28 gkg’), mi MCMHC KEE 
壤 含 盐 量 小 于 初始 值 (29.54 g-kg 和 45.79 gkg’), 
可 能 原因 是 在 特征 时 间 a,MC 和 HC 处理 2~8 cm E 
壤 含 水 量 相 比 初始 值 减少 ,并 且 2~8 cm 土 层 盐分 均 
小 于 表层 土壤 ,表明 2~8 cm 土 层 盐 分 随 着 水 分 向 上 


时 间 a~e 过 程 中 莹 发 速率 较 高 所 ,下 层 盐 分 向 上 运 
动 并 逐渐 形成 盐 结 皮 ,在 特征 时 间 d 之 后 ,蒸发 速率 
明显 减 小 ,此 时 水 分 运动 主要 以 水 汽 扩散 的 形式 向 
上 运动 ,表层 含 盐 量 趋 于 稳定 ,而 MC 和 HC 表层 土 
壤 较 为 湿润 ,水 分 主要 以 液态 水 的 形式 传输 ,携带 
盐分 向 上 缓慢 迁移 并 结晶 析出 ,表层 含 盐 量 持续 稳 
定 减 小 ,试验 结束 时 ,MC 和 HC 表层 土壤 含 盐 量 相 
比 初始 值 分 别 减 小 了 30.7% 和 15.9%。 

对 于 上 层 土壤 , 含 盐 处 理 盐 分 在 该 层 变化 幅度 
较 大 , 且 土 壤 剖 面 最 小 含 盐 量 位 于 该 层 。 在 特征 时 
HJ a~c 阶 段 ,LC 最 小 含 盐 量 所 处 深度 在 2~6 cem, MC 
和 HC 最 小 含 盐 量 在 6~10 cm。 在 特征 时 间 d~e 阶 段 
不 同 处 理 最 小 含 盐 量 所 处 深度 趋 于 稳定 。 在 试验 
结束 时 , LC 最 小 含 盐 量 为 0.31 g'kg', 位 于 2~4 em 
深度 , 相 比 初始 含 盐 量 减 小 了 90.5% ,该 处 土壤 含 盐 
量 的 突变 与 含水 量 的 突变 相 吻 合 ,进一步 验证 了 蒸 
发 锋 移动 至 2~4 cm 深度 ,此 时 盐分 主要 在 蒸发 锋 上 
方 累积 。 而 MC 和 HC 最 小 含 盐 量 在 6~8 cm 深度 ,分 
别 为 15.87 g-kg 131.08 g:kg', 相 比 初始 含 盐 量 分 
别 减 小 了 46.3% 和 32.1%。 结 果 表 明 , 盐 结 皮 越 厚 ， 

上 层 土壤 盐分 随时 间 的 减 小 幅度 越 小 。 


迁移 至 土壤 表层 ,由 于 MC 和 HC 较 高 的 盐分 浓度 其 
表层 土壤 很 快 达到 饱和 状态 ,盐分 结晶 析出 形成 盐 
结 皮 ,导致 表层 盐分 低 于 初始 值 ,而 LC 盐 结 皮 形 成 
较 慢 , 较 低 的 盐分 浓度 使 得 下 层 盐 分 持续 向 上 迁 
移 ,导致 表层 盐分 高 于 初始 值 。 随 着 试验 的 进行 ， 
LC 表层 土壤 含 盐 量 的 变化 表现 为 先 减 小 之 后 趋 于 
稳定 ,而 MC 和 HC 表层 土壤 售 盐 量 随 着 时 间 的 推移 
逐渐 减 小 。 造 成 这 种 差异 的 原因 可 能 是 LC 在 特征 


对 于 下 层 土 壤 , 土壤 含 盐 量 随 着 深度 的 增加 逐 
渐 增 加 , 盐 结 皮 越 厚 ,土壤 含 盐 量 减 小 幅度 越 小 。 
MC 和 HC 下 层 土壤 含 盐 量 相 比 LC 具 有 显著 性 差异 
(P«0.05) ,LC 盐分 峰值 位 于 表层 土壤 ,而 MC 和 HC 
盐分 峰值 位 于 下 层 土壤 。 盐 分 分 布 产 生 较 大 差异 
主要 原因 是 LC 初 始 浓度 较 低 ,土壤 表层 盐分 达到 人 饱 
和 需要 下 层 盐分 的 补给 , 盐 结 皮 形 成 较为 缓慢 , 对 
土壤 水 盐 抑 制作 用 较 小 ,下 层 盐 分 可 持续 向 上 迁 
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图 6 不 同 处 理 土壤 含水 量 实测 值 和 模拟 值 对 比 


Fig. 6 Comparison of measured and simulated values of soil water content in different treatments 


移 ,在 试验 结束 时 ,下 层 土壤 含 盐 量 相 比 初始 值 减 
小 了 67.3% 以 上 。 而 MC 和 HC 由 于 较 高 的 初始 浓 
度 , 下 层 少 量 的 盐分 向 上 迁移 使 上 层 盐 分 很 快 达到 
饱和 状态 并 结晶 析出 E s He AR DUE ROT REL a 
整个 土壤 表面 ,进而 抑制 了 水 盐 的 向 上 传输 ,使 下 
层 土壤 含 盐 量 变化 幅度 较 小 ,试验 结束 时 MC 下 层 
25-30 cm 和 HC 下 层 20~30 em 土壤 含 盐 量 相 比 初始 
值 仅 有 很 小 的 变化 。 

盐分 均衡 可 以 反映 一 定 区 域 盐分 输入 量 和 输 
出 量 之 间 的 关系 。 本 研究 中 含 盐 处 理 输入 盐分 为 
初始 饱和 不 同 浓度 盐 溶液 ,盐分 结晶 析出 形成 的 盐 
结 皮 为 输出 盐分 。 在 试验 结束 时 ,LC 表 层 、 上 层 和 
下 层 土壤 含 盐 量 分 别 为 1.04 g- kg :0.44 g kg RI 


上 层 和 下 层 土壤 含 盐 量 分 别 为 2.13 g*kg ',64.00 g- 
kg ' f 100.93 g .kg ,输出 盐分 为 33.04 g-kg ^. ZR 
表明 ,初始 盐分 浓度 越 高 ,输出 盐分 质量 越 大 , 由 于 
不 同 含 盐 处 理 蒸发 速率 的 差异 ,初始 浓度 越 高 蒸发 
速率 越 小 ,下 层 盐分 同 表层 迁移 速率 越 小 ,导致 
输出 的 盐分 占 总 盐分 质量 的 比例 越 小 。 
2.3 不 同 处 理 土壤 剖面 水 盐 动态 变化 模拟 

通过 HYDRUS-1D 模型 对 不 同 处 理 含水 量 随 时 
间 的 变化 进行 了 模拟 ,模拟 结果 如 图 6, 模 拟 精度 见 
表 2。 由 表 2 可知,CK 在 特征 时 间 b 土 壤 售 水 量 的 
模拟 精度 最 高 , 忆 和 RMSE 分 别 为 0.99 和 0.006 em 
cm”。LC 在 特征 时 间 c 土壤 含水 量 模拟 精度 最 高 ， 
尼 和 RMSE 分 别 为 0.99 和 0.006 cm-em?, MC II HC 


2.13 g* ke ', 输 出 盐分 为 7.26 g- kg; MC RE, EJ 
和 下 层 土壤 含 盐 量 分 别 为 5.30 g* kg '、17.86 g* kg ' 
和 33.20 g.kg ,输出 盐分 为 24.36 g.kg';HC 表 层 、 


土壤 含水 量 模拟 值 随 着 时 间 的 推移 模拟 精度 逐渐 
变 高 ,在 特征 时 间 e, 尼 分 别 为 0.98 和 0.99,RMSE 均 
为 0.008 cm.cm3。 从 图 6 可 以 看 出 随 着 初始 浓度 
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表 2 不 同 处 理 水 分 模拟 精度 评价 
Tab. 2 Evaluation of moisture simulation accuracy of different treatments 
不 同 处 理 R RMSE/cm'* em? 
a b c d e a b c d e 
CK 0.96 0.99 0.93 0.93 0.89 0.014 0.006 0.013 0.014 0.017 
LC 0.92 0.97 0.99 0.95 0.92 0.021 0.010 0.006 0.009 0.011 
MC 0.94 0.98 0.90 0.97 0.98 0.006 0.007 0.017 0.009 0.008 
HC 0.79 0.86 0.97 0.98 0.99 0.004 0.010 0.010 0.011 0.008 
HE ROUGE REL RMSE 为 均 方 根 误差 ;a~e 为 不 同 特征 时 间 。 下 同 。 
RI 含 盐 处 理 不 同时 间 溶 质 运 移 参数 及 模拟 精度 
Tab.3 Solute transport parameters and simulation accuracy at different times of salt treatments 
— Diem?" R RMSE/g* kg" 
0~a 阶 段 a~c 阶段 c~e 阶段 0~a 阶 段 a~c 阶 段 c~e 阶段 0~a 阶 段 a~c 阶段 c~e 阶 段 
LC 5.0 0.5 0.5 0.98 0.98 0.97 0.081 0.063 0.048 
MC 22.0 15.0 25.0 0.72 0.96 0.94 0.440 0.128 0.369 
HC 30.0 25.0 25.0 0.55 0.85 0.90 0.102 0.720 0.304 


注 :D 为 弥散 系数 。 


的 增加 ,HYDRUS-1D 模 型 对 试验 初期 土壤 含水 量 
模拟 精度 较 小 。 土 壤 含水 量 的 动态 变化 采用 VG 模 
型 进行 模拟 ,通过 分 析 V C 模型 参数 发 现 ,不 同 初始 
浓度 对 参数 a 和 nn 的 影响 较 大 。 随 着 初始 浓度 的 增 
加 ,a 逐渐 减 小 ,n 值 逐渐 增 大 ,而 a 和 nn 是 与 土壤 进 
气 值 有 关 的 参数 光 , 根 据 公式 (2) 计 算出 CK LC, 
MC fll HC 土壤 进 气 值 分 别 为 16.66.42.82 .157.97 和 
170.45, 可 以 看 出 随 着 初始 浓度 的 增加 ,土壤 进 气 值 
增 大 ,土壤 中 的 水 分 较 难 排出 ,土壤 含水 量变 化 越 
小 ,这 可 能 是 导致 试验 初期 含水 量 模拟 精度 较 差 的 

HYDRUS-1D 模型 盐分 模块 没 能 很 好 地 反映 盐 
分 在 表面 形成 结 皮 之 后 的 盐分 迁移 状况 ,因此 采用 
分 段 模拟 的 方法 ,以 更 好 模拟 盐 结 皮 土 壤 盐分 动态 
变化 过 程 , 根 据 盐 结 皮 在 不 同 蒸发 阶段 的 形成 过 程 
及 作用 机 制 纪 将 试验 过 程 分 为 3 个 阶段 ,0~a 阶 段 盐 
结 皮 开 始 形 成 并 和 覆盖 土壤 表面 ,a~c 阶段 盐 结 皮 逐 
渐 生 长 变 厚 阶 段 ,c~e 阶 段 盐 结 皮 稳 定 阶段 。 图 7 所 
示 为 含 盐 处 理 土壤 含 盐 量 的 模拟 值 ,模拟 精度 和 人 参 
数 见 表 3, 其 中 LC 各 阶段 以 及 MC 和 HC a~e 阶段 、 
c~e 阶 段 模拟 精度 较 好 , RISK 0.85, Tfj MC AI HC 
0~a 阶 段 模拟 精度 相对 较 差 。 对 比 含 盐 处 理 弥散 系 
BD) , 随 着 初始 盐分 浓度 的 增加 , D 逐渐 增 大 ,可 
能 原因 是 不 同 盐分 处 理 其 溶液 流速 不 同 ,从 而 导致 
D 产 生 较 大 差异 。 对 比 含 盐 人 处理 不 同 阶 段 ,在 试验 


a~c、c~e 阶 段 D 小 于 0~a 阶 段 , 可 能 是 试验 后 期 盐 结 
皮 改 变 了 土壤 筷 际 度 ,从 而 对 溶液 D 产 生 影响 。 

由 于 试验 条 件 限制 ,试验 仅 观测 到 MC 和 HC 处 
理 0~2 cm 土 层 的 含水 量 ,无 法 明确 蒸发 锋 是 否 下 
移 。 因 此 ,通过 HYDRUS-1D 模 型 对 表层 0.1 em T 
度 土 壤 仿 水量 进行 模拟 ,在 试验 结束 时 MC 和 HC 土 
壤 含 水 量 分 别 为 0.143 cm cem “和 0.185 em? * em ,2 
^ Ath ETE 0.1 em 处 土壤 含水 量 远大 于 风干 土 含水 
量 。 结 果 表 明 ,在 试验 周期 内 MC 和 HC 土壤 在 
0.1 em 及 以 下 深度 仍 保持 着 水 力 连 接 , 未 发 现 蒸发 
锋 的 明显 下 移 。 

通过 模型 延长 蒸发 时 间 ,继续 对 MC 和 HoC 土壤 
含水 量 进 行 监 测 。 结 果 发 现 ,MC 在 71 d 时 2 cm 深 
度 土 壤 含水 量 发 生 突 变 , 此 时 蒸发 锋 位 于 该 深度 ， 
并 且 随 独 时 间 的 推移 燕 发 锋 继 续 下 移 ,在 80 d INE 
AE PRB 4 cm 深度 ,90 d FREI 6 cm 深度 ,在 
100 d Z&  £E SEAS fa E FE 15 cm 深度 。 对 于 HC, 模 
拟 时 间 到 100 d 时 ,发 现 蒸 发 锋 位 于 土壤 2 cm 深度 ， 
120 d 蒸 发 锋 下 移 并 稳定 在 8 cm 深度。 上述 结 果 表 
明 随 着 盐 结 皮 厚 度 的 增加 ,表层 土壤 变 干 的 速度 明 
显 减缓 ,蒸发 锋 下 移 的 时 间 被 极 大 的 延长 。 

为 明确 含 盐 处 理 蒸发 锋 区 域 土 壤 含水 量 的 具 
体 变 化 ,通过 模型 对 含 盐 处 理 土 壤 含水 量 发 后 突变 
的 土 层 按 0.1 cem 划 分 剖面 网 格 对 其 土壤 含水 量 进行 
模拟 ,发 现 LC、MC 和 HoC 土壤 含水 量 突变 的 土 层 区 
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= 0~a 阶段 实测 值 —— 0~a 阶段 模拟 值 
a arc 阶段 实测 值 — a-c 阶段 模拟 值 
e c~e 阶段 实测 值 — c~e 阶段 模拟 值 


a Ova 阶段 实测 值 一 0~a 阶段 模拟 值 
^ a-c 阶段 实测 值 —— a-c 阶段 模拟 值 
e cre 阶段 实测 值 — c~e 阶段 模拟 值 


= 0~a 阶段 实测 值 — 0~a 阶段 模拟 值 
a arc 阶段 实测 值 —— a-c 阶段 模拟 值 
+ cre 阶段 实测 值 — c~e 阶段 模拟 值 


图 7 含 盐 处 理 土壤 含 盐 量 模拟 值 与 实测 值 对 比 


Fig. 7 Comparison of simulated and measured soil salt content in saline treatments 


域 仅 为 1 cm 左右 。 结 采 表 明 , BIER S pc BETETE 
差异 ,但 含 盐 处 理 燕 发 锋 区 域 均 较 罕 (1~2 em) ,该 
结果 与 Konukcu 等 "提出 的 较 宽 蒸 发 过 渡 区 的 观点 
不 一 致 ,试验 与 模型 的 结果 均 验 证 了 Gran 等 ”提出 
的 较 罕 格 发 锋 的 观点 。 


3 结论 


本 研究 通过 室内 模拟 试验 分 析 了 不 同 盐 结 皮 
厚度 对 土壤 水 盐分 布 特征 的 影响 ,并 采用 HY- 
DRUS-1D 对 土壤 水 盐 动 态 变化 过 程 进行 了 模拟 。 
研究 结果 表明 , 盐 结 皮 抑 制 了 土壤 水 盐 的 向 上 传 
输 ,并 且 盐 结 皮 越 厚 抑 制作 用 越 强 。 

(1) 对 比 无 盐 处 理 , 含 盐 处 理 盐 结 皮 越 厚 ,土壤 
剖面 含水 量 越 大 ,试验 结束 时 ,4.5 mm 厚度 盐 结 皮 
土壤 含水 量 分 布 特征 和 无 盐 处 理 相 似 ,6.6 mm 和 
7.3 mm 厚度 盐 结 皮 土壤 训 面 含水 量 显著 大 于 无 盐 
处 理 (P<0.05)。 

(2) 含 盐 处 理 随 着 盐 结 皮 厚 度 的 增加 土壤 含 
盐 量变 化 幅度 越 小 ,试验 结束 时 4.5 mm 6.6 mm 和 
7.3 mm 厚度 盐 结 皮 土 壤 放 面 最 小 含 盐 量 相 比 初始 
含 盐 量 分 别 减 小 了 90.5% 46.3% Fil 32.1%. 

(3) 通过 试验 与 模型 联合 分 析 发 现 盐 结 皮 越 
厚 ,蒸发 锋 下 移 所 需 时 间 越 长 。 本 研究 验证 了 盐 结 
皮 厚 度 会 对 土壤 水 盐分 布 特征 产生 较 大 影响 ,因此 


建议 未 来 对 于 盐 结 皮 土 壤 水 盐分 布 的 定量 分 析 综 
合 考虑 盐 结 皮 厚度 的 影 啊 。 
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Effect of salt crust thickness on distribution characteristics 
of soil water and salt 
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Abstract: Soil salt crust has an important impact on soil evolution and ecohydrological processes in arid areas. 
There are few recent studies on water and salt distribution characteristics in salt-crust soils, and the influence of 
salt-crust thickness is not considered, leading to great differences in research results. Therefore, in this paper, four 
initial salt concentration treatments (0 g-L ', 10 g*L ', 150 g-L ', and 250 g- L `) were set to obtain different salt- 
crust thicknesses (4.5 mm, 6.6 mm, 7.3 mm) through laboratory simulation tests, and the soil-profile dynamics of 
water and salt were compared and analyzed using a partial repeated stepwise withdrawal method. The results 
were as follows: (1) Compared with the non-salt treatment, the thicker the salt crust, the larger the soil-profile wa- 
ter content, and the smaller the salt- content variation range. (2) At the end of the experiment, the water content 
distribution characteristics of the 4.5 mm salt-crust soil were similar to those of the unsalted treatment, and the 
water contents of the 6.6 mm and 7.3 mm salt-crust soils were significantly higher than that of the unsalted treat- 
ment (P«0.05). (3) At the end of the test, the minimum salt contents of the 4.5 mm, 6.6 mm, and 7.3 mm salt- 
crust soils decreased by 90.5%, 46.3%, and 32.1%, respectively, compared with their initial salt contents. The re- 
sults confirm that salt-crust thickness has a great influence on the distribution of soil water and salt. Therefore, it 
is suggested that the influence of salt-crust thickness should be considered comprehensively in future research on 
distribution characteristics of water and salt. 

Key words: salt crust thickness; distribution characteristics of water and salt; evaporation front; HYDRUS-1D 


model 


